Numerical Simulation of the Nernst Effect in Extreme Type-II Superconductors: 
A Negative Nernst Signal and its Noise Power Spectra 



(N 

O 
(N 

o 

Q 

oo 



Sangwoo S. Chung/' Paata Kakashvili, 

3.1.0 and C. J. Bolech^'i 

^Physics & Astronomy Department, Rice University, 6100 Main Street, Houston, TX 71005 
^Department of Physics, University of Cincinnati, 345 Clifton Court, Cincinnati, OH 45221 
^NORDITA, Roslagstullsbacken 23, 106 91 Stockholm, Sweden 
(Dated: January 1, 2013) 

Recently, different transport coefficients have been measured in liigli- Tc superconductors to pin- 
point tlie nature of tfie pseudogap pfiase. In particular, the thermoelectric coefficients received 
considerable attention both theoretically and experimentally. We numerically simulate the Nernst 
effect in extreme type-II superconductors using the time-dependent Ginzburg-Landau equations. 
We report the sign reversal of the thermoelectric coefficient axy at temperatures close to the mean- 
field transition temperature (H), which qualitatively agrees with recent experiments on high- Tc 
materials. We also discuss the noise power spectrum of Oxy, which shows 1//'' behavior. Based 
on this observation, we propose an experiment to distinguish among different regimes of vortex 
dynamics by measuring the noise correlations of the Nernst signal. 
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I. INTRODUCTION 



The unexpected observation of unusually large Nernst 
coefficients (v) in cuprates at temperatures much greater 
than have drawn much spotlight to the Nernst ef- 
fect over the past decade. As the Nernst coefHcients are 
generally small in most ordinary metals,'^ the large val- 
ues of v seen at T > Tc have elevated the Nernst ef- 
fect as one of the central avenues used in attempts to 
characterize the nature of the pseudogap state. Since 
then, an extensive investigation on the subject has been 
done, both experimentallji^"— and theoretically,—"— 
producing widely differing proposals on the origin of 
the phenomenon. Most of these competing interpreta- 
tions center in the dynamics of either vorticesi^"— i^i"— 
or quasiparticlesji"— 

In this article, we use the time-dependent Ginzburg- 
Landau (TDGL) equations to simulate the Nernst effect 
for an extreme type-II superconductor and (i) report a 
sign reversal of the transverse thermoelectric coefficient 
axy consistent with the one observed in experiments, and 
(ii) examine the noise properties of the Nernst signal, 
which has been surprisingly overlooked thus far for the 
Nernst effect and, we find, can distinguish among differ- 
ent dynamical phases of the vortex system. 

The Nernst effect is the emergence of a transverse volt- 
age V under the presence of an applied thermal gradi- 
ent Vr||y and a perpendicular magnetic field H||z (see 
Fig. [1]). From linear- response theory, the transport elec- 
tric current due to an applied electric field E and a ther- 
mal gradient VT has the form J = (tE — dVT, where 
a and a are the electric and thermoelectric conductiv- 
ity tensors. Restricted to systems with negligible Hall 
effect, the off-diagonal term of a is related to the Nernst 
coefficient as v{H,T) w axy{H,T)/ Haxx{H,T)^ 




FIG. 1. (Color online) Geometry of a Nernst experiment. 
A transverse voltage V in the x direction appears when a 
magnetic field H and a thermal gradient VT are applied si- 
multaneously in the z and y directions, respectively. 



II. EXTREME TYPE-II LIMIT 

The TDGL equations can formally be derived by tak- 
ing the Ginzburg-Landau free-energy functional2£ and 
employing Langevin dynamics for the time evolution of 
the superconducting order parameter ip. For quintessen- 
tial high-Tc superconductors, it is appropriate to con- 
sider the extreme type-II hmit where k, the ratio of the 
penetration depth A and the coherence length ^, is as- 
sumed to be infinite by taking A — ?> oo. In this sim- 
pUfied picture, the applied magnetic field in the sam- 
ple is uniform and unscreened, as the lower critical field 
Hci ^ {Iiik)/k 022. Therefore, H(r,t) = ffext(Oz and 
we select an instantaneous Landau gauge for the vector 
potential A(r,i) = iJext(i) (r x z) /2. Consequently, the 
V X H term in Maxwell's equations vanishes, and we ob- 
tain the following relations that govern the evolution of 
■(/'(r, t) and the electric potential 0(r, t), 



{dt + icj)) ip 



dtA + V0 : 



--\(-iV - Afi; 

7] 

+(i-r)(|V'p-i)^]+/, (1) 

(1 -r)Re[V'*(-iV- A)i/;], (2) 
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along with zero perpendicular supercurrent boundary 
conditions, (— iV — A)-0|n — 0. 

Temperature, length and time are measured in 
units of the mean- field transition temperature , 
zero-temperature coherence length ^(0), and to = 
TTh/{96kBT^^), respectively. Vector potential A and 
magnetic field H are in units of i?c2 (0)^(0) and Hc2{0), 
respectively, where i?c2(0) is the upper-critical field at 
T ~ 0. The dimensionless parameter ij is proportional to 
the ratio of characteristic times for "0 and A, and / is a 
random thermal noise. The order parameter ^ has been 
normalized by the equilibrium value ipQo{T). 



III. NUMERICAL SCHEME 

Initially inspired by the Euler-forward based finite- 
difference scheme in Refs. |21| and i22 (cf. Ref. i23l ). 
we develop a gauge-invariant discrete implementation 
of Eqs. ([Ij and ©. Space and time have been 
discretized as shown in Fig. [5] with ipij^nS placed 
on the nodes, which are at fixed local temperatures 
Tij. In addition to the spacelike link variables. 



dS^Ay{xi,^,tn)), we introduce a time-like 



link variable Wi 



exp 



corporate (f> into the numerics. 
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FIG. 2. (Color online) Scheme of discretization. The or- 
der parameter i/i is a nodal variable and , , and W are 
placed on links connecting the nodes. We take the grid spac- 
ing Ax = Ay = a = 0.5, while we set At — 0.015 for the time 
interval between successive time steps, both in normalized 
units. These values were chosen so as to be able to resolve 
the structure and dynamics of the vortices. Use of smaller 
values does not alter the results in any significant way. 

To write Eqs. (dJ and © in discrete forms, accurate 
to order O(Ai^) and O(a^), we first note that 



exp i / dr [drA^ + d^cj)] 



v=yj 



1 + iaAt [dtA.^ + 9^ 



where Xi+1/2 = {xi + Xt+i)/2. Using Eq. ^ 
and the relation Reltp* {—idx — Ax)i']^,_^^^^ y_ « 

ilm (V'*,C/i^^-^,+i,,)^, Eq. m yields 



l + iAt{l 



(4) 



and together with an analogous expression in the y 
direction, they provide a solution to a Neumann-type 



[with ■(/'i,j,n's given. 



m 



boundary-value problem for Wij^ 
and J7ij,„'s fixed by ^^oxt(tn)]- 

Moreover, it can be shown that Eq. (HI) — 
terms of the discrete variables can be expressed as 

(Wij^j^nViJ.n+l - / At = Vij^n + /ij,„ with 

-(l-T,,,)(V^,,V'-j-l)^^,]„, (5) 

from which ipij's at t = tn+i may be evaluated. 

Wc perform our simulations on square grids with N xN 
cells with N ranging from 64 to 384. A time-independent 
thermal gradient is realized by linearly varying the lo- 
cal temperature Ti,j along the y direction of the sam- 
ple. For N = 192, we apply a fixed temperature dif- 
ference AT = 0.03 throughout, so that temperatures at 
the coldest and the hottest edges are Tip = T =F AT/2, 
respectively, where T marks the average temperature of 
the sample. The thermal fluctuation term / is a Gaus- 
sian noise with standard deviation aij — iiEoTi^jAt/Q, 
where Eg is a material-dependent parameter. For our 
choice of 77 = 1, we set Eq = \ (guided by the fluctuation- 
dissipation theorem and values used in the literature to 
fit experiments-^^) and confirmed that the results do not 
depend sensitively on variations around this value. 



IV. SIMULATION RESULTS 



Our initial setting at i = is the Meissner state, where 
V'ij = Uf^j = U^j — 1. With the fixed thermal gradient 
applied, we ramp up the magnetic field to a set value 
in the first few thousand steps and wait for 0.2 million 
time steps to reach a steady state. We then time average 
ctxy{tn) over 0.8 million time steps to obtain Uxy, where 



Jx{tn) _ J J Jx{x,y,tn)dxdy 



IVTI 



NaAT 



(6) 



and Jx^{l- T)Re [ip*{~idx ~ Ax)^^] at t = t^. 



A. Transverse thermoelectric response 



i + l/2,yj,*n ' 



Figure [3l[a) is the plot of axy vs temperature for var- 
ious magnetic fields (see Fig. [S] for axy vs H). At low 
(3) temperatures and, particularly, at lower magnetic fields. 
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axy attains large positive values as also seen in both 
experimental^ and theoretica l^^i^^ results. As T is in- 
creased, axy falls monotonically until it reaches a vortex- 
proliferation temperature Typ{H), which can be iden- 
tified from Fig. ^h) (see Fig. [5] for vortex densities vs 
H). At Tvp{H), the number of vortices and antivortices 
proliferate and the long-range phase coherence of ip is de- 
stroyed. The exhibition of high axy values at T < Typ 
is consistent with the very low resistances of the samples 
in that regime (cf. Ref. [23). Above Typ{H), the curves 
of axy resemble the profiles seen in Nernst coefficients 
for typical high- Tc superconductors. At even higher tem- 
peratures, axy becomes negative, and we designate the 
onset temperature for negative Nernst signals as Tnn (H) 
for later discussions. 

Above Tnn, we find that the onset of negative Nernst 
signal is correlated to the onset of the superconducting 
to normal transition at the hotter edge. Prior to the 
transition, H lies in the range Hc2{T-) < H < Hc3{T+), 
where iJcS is the critical field that destroys surface su- 
perconductivity. The bulk is already at the normal state, 
while the Nernst signal is dominated by superconducting 
edge currents circulating around the bulk of the sam- 
ple. As T is continually increased, Hcs (T) decreases and 
the hotter edge crosses the phase boundary to enter the 
normal state. Consequently, the Cooper-pair density at 
that edge is abruptly suppressed, resulting in a signifi- 
cant asymmetry in the supercurrent density relative to 
the colder edge and a different vortex vs antivortex dy- 
namics across the sample. 

From experiments, the negative axy reported in Ref.fiol 
for an optimally doped YBCO sample displays a notable 
resemblance to our results. A qualitatively similar profile 
of negative Nernst signals can also be inferred from mea- 
surements of Nernst coefficients done on a wide range 
of samples at various doping levels,*'^ which indicates 
that this effect is ubiquitous and independent of details 
or fine tuning. Theoretical results on granular super- 
conductors in Ref. [l^ also show negative Nernst signals, 
which they interpret as diffusion of vortex vacancies in 
a pinned-vortex state. In our simulations, however, the 
sign reversal is not from quasiparticle effects as that 
scenario is not comprised in the Ginzburg-Landau frame- 
work. Moreover, vortex vacancies do not seem to be the 
cause, since it can be inferred from Ref. [Hthat Tnn (if) 
would increase as H is increased, while in Fig. [3Ka) we 
see that Tt^^{H) displays an opposite trend as the field 
is increased. 



B. Noise analysis 

The dynamic phase diagram of vortex matter is more 
diverse than the one at equilibrium.-^ In our simula- 
tions, at finite k and in the absence of pinning centers, 
we observe rich dynamical behavior. Depending on the 
strength of the driving force (temperature gradient in 
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FIG. 3. (Color online) (a) Temperature dependence of a^y 
Inset: a^y zoomed into the negative-Nernst-signal region with 
TNN(Jy)'s marked by vertical arrows; the error bars (com- 
parable to symbol sizes) highlight the statistical accuracy of 
the negative signals. The error bar for each data point has 
been obtained from eight independent measurements of Qa;y, 
each measurement consisting of time averaging over 0.1 mil- 
lion time steps, (b) Temperature dependence of n^+v, the 
total density of vortices plus antivortices, showing the onset 
of the vortex-prohferation regime, which is correlated to small 
kinks in axy and is seen to be unrelated to the onset of the 
negative Nernst signal. The net vorticity n^-v (dotted lines) 
reflects the applied magnetic flux that is independent of tem- 
perature. 



our case), temperature, and magnetic field, vortex mo- 
tion can be characterized into different regimes: immo- 
bile lattice of vortices due to the geometrical pinning in 
the finite-size sample, sliding vortex crystal, and flowing 
vortex liquids. It is not surprising that noise properties 
of transport signals depend on the dynamical behavior of 
vortices. Especially, since the Nernst signal is highly sen- 
sitive to vortex dynamics, we propose the measurement 
of its noise properties as a means to tell apart different 
dynamical scenarios. For this purpose, we have computed 
the noise correlation for axy{t). The noise power S{f) of 
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a signal x{t) is defined as the Fourier transform of the 
temporal correlation function, 

Sit (Ax{t)Ax{t') + Ax{t')Ax{t)) , (7) 

where Ax{t) = x{t) — {x). 

Depending on the dynamical regimes described above, 
S{f) shows different qualitative behaviors. We note here 
that in the k — > oo limit, the shear modulus of the vor- 
tex matter vanishes and thus only the vortex liquid phase 
survives. To give a full picture, we provide a short de- 
scription of the noise behavior in the finite- k regime also, 
but a more detailed account will be published elsewhere. 
In the pinned phase, the Nernst signal is suppressed and 
the noise power is nearly independent of /, thus display- 
ing a white-noise spectrum. In the sliding-lattice phase, 
the noise power shows a characteristic narrow-band be- 
havior and most of the spectral density is found near a 
particular frequency. The characteristic frequency can be 
estimated as ^ v/£, where v is the speed of vortices 
and £ is the distance between vortices in the direction of 
motion. This type of noise is known as washboard noise 
and has been observed experimentally in the conduction 
noise spectra of driven BSCCO samples. In the vortex- 
liquid regime, noise power again shows a wide-band be- 
havior. In the infinite-K limit, we found it to be fitted 
well by a power law 1 //^ . Figure |4] shows the plot of (3 
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FIG. 4. (Color online) Exponent of the noise power, j3, as a 
function of T for different values of H. The error bars (shown 
for selected data points) have been obtained from a least- 
squares fitting procedure on the log S vs log / plots. The 
vertical arrows mark the onset temperatures (Tnn) for the 
negative Nernst signal obtained using linear interpolation in 
Fig. 3(a). 

for different magnetic fields as a function of temperature. 
We observe that in the low-T vortex-liquid regime, the 
exponent /3 « 1 and is nearly independent of the applied 
magnetic field. As the system enters into the negative- 
Nernst-signal regime, /3 drastically drops to lower values. 
This sharp feature provides an excellent means for ex- 



periments to causally correlate the onset temperature for 
negative Nernst signals, Tnn, to vortex dynamics. As the 
temperature approaches the mean-field transition tem- 
perature, the exponent is further reduced, but remains 
finite. In contrast to the low-T vortex-liquid regime, in 
this case, (3 strongly depends on the applied magnetic 
field. We remark that we do not see any observable signa- 
tures of the vortex-proliferation transition in the behavior 
of /3. Even though 1// behavior of noise is ubiquitous in 
transport measurements of many physical systems to 
the best of our knowledge, it has not been measured for 
the Nernst signal. In YBCO samples, a similar power-law 
behavior has been observed in flux noise spectra in a non- 
driven system2S and in longitudinal- voltage noise spectra 
in a current-driven system. "^^ Although the above mea- 
surements have been done for regimes and observables 
other than those presented here, they illustrate the feasi- 
bility of the noise-measurement scenario we propose here. 



V. CLOSING REMARKS 

In this article, a computational method has been in- 
troduced to simulate the Nernst effect in high- Tc su- 
perconductors using the TDGL theoretical framework. 
We show that the thermoelectric coefficient axy (i) is 
large and monotonically decreasing for T < Typ , (ii) dis- 
plays a profile that is qualitatively consistent with ex- 
periments at intermediate T's, and (iii) reverses signs 
at r > Tnn > 2vp- Furthermore, we propose a noise- 
based observable to distinguish between different dynam- 
ical phases of the vortex system. The noise power of the 
Nernst signal reveals distinctive qualitative and quantita- 
tive properties for different regimes, thus demonstrating 
its great potential as an experimental probe of vortex 
dynamics. 
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APPENDIX 

For ease of comparison with the existing literature, we 
alternatively present the data in Figs. [Sja) and ^h) as 
functions of magnetic fields for several temperatures [con- 
trary to functions of temperatures for several fields as in 
Figs, ^a) and[3]Jb)]. For any given curve in Fig. [SI we 
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see a linear increase of vortex densities until reaching a 
certain field value, at which vortices proliferate due to 
the onset of vortex-antivortex pair generation. The pro- 
liferation starts at higher fields for lower temperatures. 



from the negative side. The sign reversal of Uxy occurs 
at higher magnetic fields and the magnetic-field ranges 
for negative axy broaden for lower temperatures. 
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FIG. 5. (Color online) Magnetic field dependence of n^+v, the 
total density of vortices plus antivortices, showing the onset 
of the vortex-proliferation regime. 

Figure [S] shows the plot of axy in the region where 
the negative values are conspicuous. As H is increased 
beyond Hc2{T) for a given temperature, axy reverses sign 
and becomes negative before diminishing towards zero 
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FIG. 6. (Color online) Magnetic field dependence of Oxy', 
the error bars (comparable to symbol sizes) highlight the sta- 
tistical accuracy of the negative signals. The error bar for 
each data point has been obtained from eight independent 
measurements of a^y, each measurement consisting of time 
averaging over O.f million time steps. 
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